##Add Bayesian tests functions
#create function to conduct the Bayesian Sign Test
BayesianSignTest <- function(diffVector,rope_min,rope_max) {
library(MCMCpack)
samples <- 3000
#build the vector 0.5 1 1 ....... 1
weights <- c(0.5,rep(1,length(diffVector)))
#add the fake first observation in 0
diffVector <- c (0, diffVector)
#for the moment we implement the sign test. Signedrank will follows
probLeft <- mean (diffVector < rope_min)
probRope <- mean (diffVector > rope_min & diffVector < rope_max)
probRight <- mean (diffVector > rope_max)
results = list ("probLeft"=probLeft, "probRope"=probRope,
"probRight"=probRight)
return (results)
}
##Create function to conduct Bayesian Signed Rank Test
BayesianSignedRank <- function(diffVector,rope_min,rope_max) {
library(MCMCpack)
samples <- 30000
#build the vector 0.5 1 1 ....... 1
weights <- c(0.5,rep(1,length(diffVector)))
#add the fake first observation in 0
diffVector <- c (0, diffVector)
sampledWeights <- rdirichlet(samples,weights)
winLeft <- vector(length = samples)
winRope <- vector(length = samples)
winRight <- vector(length = samples)
for (rep in 1:samples){
currentWeights <- sampledWeights[rep,]
for (i in 1:length(currentWeights)){
for (j in 1:length(currentWeights)){
product= currentWeights[i] * currentWeights[j]
if (diffVector[i]+diffVector[j] > (2*rope_max) ) {
winRight[rep] <- winRight[rep] + product
}
else if (diffVector[i]+diffVector[j] > (2*rope_min) ) {
winRope[rep] <- winRope[rep] + product
}
else {
winLeft[rep] <- winLeft[rep] + product
}
}
}
maxWins=max(winRight[rep],winRope[rep],winLeft[rep])
winners = (winRight[rep]==maxWins)*1 + (winRope[rep]==maxWins)*1 + (winLeft[rep]==maxWins)*1
winRight[rep] <- (winRight[rep]==maxWins)*1/winners
winRope[rep] <- (winRope[rep]==maxWins)*1/winners
winLeft[rep] <- (winLeft[rep]==maxWins)*1/winners
}
results = list ("winLeft"=mean(winLeft), "winRope"=mean(winRope),
"winRight"=mean(winRight) )
return (results)
}
#Create function to conduct the Bayesian Correlated t.test
#diff_a_b is a vector of differences between the two classifiers, on each fold of cross-validation.
#If you have done 10 runs of 10-folds cross-validation, you have 100 results for each classifier.
#You should have run cross-validation on the same folds for the two classifiers.
#Then diff_a_b is the difference fold-by-fold.
#rho is the correlation of the cross-validation results: 1/(number of folds)
#rope_min and rope_max are the lower and the upper bound of the rope
correlatedBayesianTtest <- function(diff_a_b,rho,rope_min,rope_max){
if (rope_max < rope_min){
stop("rope_max should be larger than rope_min")
}
delta <- mean(diff_a_b)
n <- length(diff_a_b)
df <- n-1
stdX <- sd(diff_a_b)
sp <- sd(diff_a_b)*sqrt(1/n + rho/(1-rho))
p.left <- pt((rope_min - delta)/sp, df)
p.rope <- pt((rope_max - delta)/sp, df)-p.left
results <- list('left'=p.left,'rope'=p.rope,'right'=1-p.left-p.rope)
return (results)
}
set.seed(16974)
###################################################5.40.5 ROPE Comparisons for Dissertation
##Random Forest Results
rf_dataset_av<-c(0.8589, 0.9240, 0.9797)
rf_pca.5.40.5_n1_av<-c(0.2413, 0.5473, 0.9629)
rf_pca.5.40.5_n2_av<-c(0.3800, 0.8208, 0.9184)
rf_pca.5.40.5_n3_av<-c(0.5292, 0.4373, 0.1982)
rf_pca.5.40.5_n4_av<-c(0.8174, 0.3127, 0.1082)
rf_pca.5.40.5_n5_av<-c(0.7592, NA, NA)
rf_kde.5.40.5_n1_av<-c(0.9074, 0.9135, 0.9909)
rf_kde.5.40.5_n2_av<-c(0.9068, 0.8922, 0.9348)
rf_kde.5.40.5_n3_av<-c(0.9036, 0.7078, 0.9717)
rf_kde.5.40.5_n4_av<-c(0.8727, 0.4971, 0.9609)
rf_kde.5.40.5_n5_av<-c(0.8489, 0.4191, 0.9296)
######################## ROPE PCA
diff_rf_pca.5.40.5_n1_av<-rf_dataset_av - rf_pca.5.40.5_n1_av
bsr_diff_rf_pca.5.40.5_n1_av<-BayesianSignedRank(as.matrix(diff_rf_pca.5.40.5_n1_av),-0.01,0.01)
bsr_diff_rf_pca.5.40.5_n1_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.03726667
##
## $winRight
## [1] 0.9627333
bsr_diff_rf_pca.5.40.5_n1_av_odds.left<-bsr_diff_rf_pca.5.40.5_n1_av $winLeft/bsr_diff_rf_pca.5.40.5_n1_av $winRight
bsr_diff_rf_pca.5.40.5_n1_av_odds.left
## [1] 0
plot(rope(diff_rf_pca.5.40.5_n1_av,c(-0.01,0.01)))

diff_rf_pca.5.40.5_n2_av<-rf_dataset_av - rf_pca.5.40.5_n2_av
bsr_diff_rf_pca.5.40.5_n2_av<-BayesianSignedRank(as.matrix(diff_rf_pca.5.40.5_n2_av),-0.01,0.01)
bsr_diff_rf_pca.5.40.5_n2_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.0087
##
## $winRight
## [1] 0.9913
bsr_diff_rf_pca.5.40.5_n2_av_odds.left<-bsr_diff_rf_pca.5.40.5_n2_av $winLeft/bsr_diff_rf_pca.5.40.5_n2_av $winRight
bsr_diff_rf_pca.5.40.5_n2_av_odds.left
## [1] 0
plot(rope(diff_rf_pca.5.40.5_n2_av,c(-0.01,0.01)))

diff_rf_pca.5.40.5_n3_av<-rf_dataset_av - rf_pca.5.40.5_n3_av
bsr_diff_rf_pca.5.40.5_n3_av<-BayesianSignedRank(as.matrix(diff_rf_pca.5.40.5_n3_av),-0.01,0.01)
bsr_diff_rf_pca.5.40.5_n3_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.008533333
##
## $winRight
## [1] 0.9914667
bsr_diff_rf_pca.5.40.5_n3_av_odds.left<-bsr_diff_rf_pca.5.40.5_n3_av $winLeft/bsr_diff_rf_pca.5.40.5_n3_av $winRight
bsr_diff_rf_pca.5.40.5_n3_av_odds.left
## [1] 0
plot(rope(diff_rf_pca.5.40.5_n3_av,c(-0.01,0.01)))

diff_rf_pca.5.40.5_n4_av<-rf_dataset_av - rf_pca.5.40.5_n4_av
bsr_diff_rf_pca.5.40.5_n4_av<-BayesianSignedRank(as.matrix(diff_rf_pca.5.40.5_n4_av),-0.01,0.01)
bsr_diff_rf_pca.5.40.5_n4_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.0097
##
## $winRight
## [1] 0.9903
bsr_diff_rf_pca.5.40.5_n4_av_odds.left<-bsr_diff_rf_pca.5.40.5_n4_av $winLeft/bsr_diff_rf_pca.5.40.5_n4_av $winRight
bsr_diff_rf_pca.5.40.5_n4_av_odds.left
## [1] 0
plot(rope(diff_rf_pca.5.40.5_n4_av,c(-0.01,0.01)))

#diff_rf_pca.5.40.5_n5_av<-rf_dataset_av - rf_pca.5.40.5_n5_av
#bsr_diff_rf_pca.5.40.5_n5_av<-BayesianSignedRank(as.matrix(diff_rf_pca.5.40.5_n5_av),-0.01,0.01)
#bsr_diff_rf_pca.5.40.5_n5_av
#bsr_diff_rf_pca.5.40.5_n5_av_odds.left<-bsr_diff_rf_pca.5.40.5_n5_av $winLeft/bsr_diff_rf_pca.5.40.5_n5_av $winRight
#bsr_diff_rf_pca.5.40.5_n5_av_odds.left
#plot(rope(diff_rf_pca.5.40.5_n5_av,c(-0.01,0.01)))
########################## ROPE KDE
diff_rf_kde.5.40.5_n1_av<-rf_dataset_av - rf_kde.5.40.5_n1_av
bsr_diff_rf_kde.5.40.5_n1_av<-BayesianSignedRank(as.matrix(diff_rf_kde.5.40.5_n1_av),-0.01,0.01)
bsr_diff_rf_kde.5.40.5_n1_av
## $winLeft
## [1] 0.6762667
##
## $winRope
## [1] 0.2452667
##
## $winRight
## [1] 0.07846667
bsr_diff_rf_kde.5.40.5_n1_av_odds.left<-bsr_diff_rf_kde.5.40.5_n1_av $winLeft/bsr_diff_rf_kde.5.40.5_n1_av $winRight
bsr_diff_rf_kde.5.40.5_n1_av_odds.left
## [1] 8.618522
plot(rope(diff_rf_kde.5.40.5_n1_av,c(-0.01,0.01)))

diff_rf_kde.5.40.5_n2_av<-rf_dataset_av - rf_kde.5.40.5_n2_av
bsr_diff_rf_kde.5.40.5_n2_av<-BayesianSignedRank(as.matrix(diff_rf_kde.5.40.5_n2_av),-0.01,0.01)
bsr_diff_rf_kde.5.40.5_n2_av
## $winLeft
## [1] 0.1574
##
## $winRope
## [1] 0.2224
##
## $winRight
## [1] 0.6202
bsr_diff_rf_kde.5.40.5_n2_av_odds.left<-bsr_diff_rf_kde.5.40.5_n2_av $winLeft/bsr_diff_rf_kde.5.40.5_n2_av $winRight
bsr_diff_rf_kde.5.40.5_n2_av_odds.left
## [1] 0.2537891
plot(rope(diff_rf_kde.5.40.5_n2_av,c(-0.01,0.01)))

diff_rf_kde.5.40.5_n3_av<-rf_dataset_av - rf_kde.5.40.5_n3_av
bsr_diff_rf_kde.5.40.5_n3_av<-BayesianSignedRank(as.matrix(diff_rf_kde.5.40.5_n3_av),-0.01,0.01)
bsr_diff_rf_kde.5.40.5_n3_av
## $winLeft
## [1] 0.3029
##
## $winRope
## [1] 0.1986333
##
## $winRight
## [1] 0.4984667
bsr_diff_rf_kde.5.40.5_n3_av_odds.left<-bsr_diff_rf_kde.5.40.5_n3_av $winLeft/bsr_diff_rf_kde.5.40.5_n3_av $winRight
bsr_diff_rf_kde.5.40.5_n3_av_odds.left
## [1] 0.6076635
plot(rope(diff_rf_kde.5.40.5_n3_av,c(-0.01,0.01)))

diff_rf_kde.5.40.5_n4_av<-rf_dataset_av - rf_kde.5.40.5_n4_av
bsr_diff_rf_kde.5.40.5_n4_av<-BayesianSignedRank(as.matrix(diff_rf_kde.5.40.5_n4_av),-0.01,0.01)
bsr_diff_rf_kde.5.40.5_n4_av
## $winLeft
## [1] 0.08063333
##
## $winRope
## [1] 0.2419333
##
## $winRight
## [1] 0.6774333
bsr_diff_rf_kde.5.40.5_n4_av_odds.left<-bsr_diff_rf_kde.5.40.5_n4_av $winLeft/bsr_diff_rf_kde.5.40.5_n4_av $winRight
bsr_diff_rf_kde.5.40.5_n4_av_odds.left
## [1] 0.1190277
plot(rope(diff_rf_kde.5.40.5_n4_av,c(-0.01,0.01)))

diff_rf_kde.5.40.5_n5_av<-rf_dataset_av - rf_kde.5.40.5_n5_av
bsr_diff_rf_kde.5.40.5_n5_av<-BayesianSignedRank(as.matrix(diff_rf_kde.5.40.5_n5_av),-0.01,0.01)
bsr_diff_rf_kde.5.40.5_n5_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.03673333
##
## $winRight
## [1] 0.9632667
bsr_diff_rf_kde.5.40.5_n5_av_odds.left<-bsr_diff_rf_kde.5.40.5_n5_av $winLeft/bsr_diff_rf_kde.5.40.5_n5_av $winRight
bsr_diff_rf_kde.5.40.5_n5_av_odds.left
## [1] 0
plot(rope(diff_rf_kde.5.40.5_n5_av,c(-0.01,0.01)))

################################ Support Vector Machine
##Support Vector Machine Results
svm_dataset_av<-c(0.8270, 0.9275, 0.9769)
svm_pca.5.40.5_n1_av<-c(0.3076, 0.5336, 0.8983)
svm_pca.5.40.5_n2_av<-c(0.3135, 0.8576, 0.9045)
svm_pca.5.40.5_n3_av<-c(0.7360, 0.4341, 0.3384)
svm_pca.5.40.5_n4_av<-c(0.7733, 0.3370, 0.1017)
svm_pca.5.40.5_n5_av<-c(0.7592, NA, NA)
svm_kde.5.40.5_n1_av<-c(0.8500, 0.9047, 0.8992)
svm_kde.5.40.5_n2_av<-c(0.8548, 0.8713, 0.9022)
svm_kde.5.40.5_n3_av<-c(0.8455, 0.6194, 0.8998)
svm_kde.5.40.5_n4_av<-c(0.8140, 0.3931, 0.8994)
svm_kde.5.40.5_n5_av<-c(0.8374, 0.3005, 0.8983)
######################## ROPE PCA
diff_svm_pca.5.40.5_n1_av<-svm_dataset_av - svm_pca.5.40.5_n1_av
bsr_diff_svm_pca.5.40.5_n1_av<-BayesianSignedRank(as.matrix(diff_svm_pca.5.40.5_n1_av),-0.01,0.01)
bsr_diff_svm_pca.5.40.5_n1_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.008133333
##
## $winRight
## [1] 0.9918667
bsr_diff_svm_pca.5.40.5_n1_av_odds.left<-bsr_diff_svm_pca.5.40.5_n1_av$winLeft/bsr_diff_svm_pca.5.40.5_n1_av $winRight
bsr_diff_svm_pca.5.40.5_n1_av_odds.left
## [1] 0
plot(rope(diff_svm_pca.5.40.5_n1_av,c(-0.01,0.01)))

diff_svm_pca.5.40.5_n2_av<-svm_dataset_av - svm_pca.5.40.5_n2_av
bsr_diff_svm_pca.5.40.5_n2_av<-BayesianSignedRank(as.matrix(diff_svm_pca.5.40.5_n2_av),-0.01,0.01)
bsr_diff_svm_pca.5.40.5_n2_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.008233333
##
## $winRight
## [1] 0.9917667
bsr_diff_svm_pca.5.40.5_n2_av_odds.left<-bsr_diff_svm_pca.5.40.5_n2_av$winLeft/bsr_diff_svm_pca.5.40.5_n2_av $winRight
bsr_diff_svm_pca.5.40.5_n2_av_odds.left
## [1] 0
plot(rope(diff_svm_pca.5.40.5_n2_av,c(-0.01,0.01)))

diff_svm_pca.5.40.5_n3_av<-svm_dataset_av - svm_pca.5.40.5_n3_av
bsr_diff_svm_pca.5.40.5_n3_av<-BayesianSignedRank(as.matrix(diff_svm_pca.5.40.5_n3_av),-0.01,0.01)
bsr_diff_svm_pca.5.40.5_n3_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.008933333
##
## $winRight
## [1] 0.9910667
bsr_diff_svm_pca.5.40.5_n3_av_odds.left<-bsr_diff_svm_pca.5.40.5_n3_av$winLeft/bsr_diff_svm_pca.5.40.5_n3_av $winRight
bsr_diff_svm_pca.5.40.5_n3_av_odds.left
## [1] 0
plot(rope(diff_svm_pca.5.40.5_n3_av,c(-0.01,0.01)))

diff_svm_pca.5.40.5_n4_av<-svm_dataset_av - svm_pca.5.40.5_n4_av
bsr_diff_svm_pca.5.40.5_n4_av<-BayesianSignedRank(as.matrix(diff_svm_pca.5.40.5_n4_av),-0.01,0.01)
bsr_diff_svm_pca.5.40.5_n4_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.0093
##
## $winRight
## [1] 0.9907
bsr_diff_svm_pca.5.40.5_n4_av_odds.left<-bsr_diff_svm_pca.5.40.5_n4_av$winLeft/bsr_diff_svm_pca.5.40.5_n4_av $winRight
bsr_diff_svm_pca.5.40.5_n4_av_odds.left
## [1] 0
plot(rope(diff_svm_pca.5.40.5_n4_av,c(-0.01,0.01)))

#diff_svm_pca.5.40.5_n5_av<-svm_dataset_av - svm_pca.5.40.5_n5_av
#bsr_diff_svm_pca.5.40.5_n5_av<-BayesianSignedRank(as.matrix(diff_svm_pca.5.40.5_n5_av),-0.01,0.01)
#bsr_diff_svm_pca.5.40.5_n5_av
#bsr_diff_svm_pca.5.40.5_n5_av_odds.left<-bsr_diff_svm_pca.5.40.5_n5_av$winLeft/bsr_diff_svm_pca.5.40.5_n5_av $winRight
#bsr_diff_svm_pca.5.40.5_n5_av_odds.left
#plot(rope(diff_svm_pca.5.40.5_n5_av,c(-0.01,0.01)))
########################## ROPE KDE
diff_svm_kde.5.40.5_n1_av<-svm_dataset_av - svm_kde.5.40.5_n1_av
bsr_diff_svm_kde.5.40.5_n1_av<-BayesianSignedRank(as.matrix(diff_svm_kde.5.40.5_n1_av),-0.01,0.01)
bsr_diff_svm_kde.5.40.5_n1_av
## $winLeft
## [1] 0.1544667
##
## $winRope
## [1] 0.04713333
##
## $winRight
## [1] 0.7984
bsr_diff_svm_kde.5.40.5_n1_av_odds.left<-bsr_diff_svm_kde.5.40.5_n1_av $winLeft/bsr_diff_svm_kde.5.40.5_n1_av $winRight
bsr_diff_svm_kde.5.40.5_n1_av_odds.left
## [1] 0.1934703
plot(rope(diff_svm_kde.5.40.5_n1_av,c(-0.01,0.01)))

diff_svm_kde.5.40.5_n2_av<-svm_dataset_av - svm_kde.5.40.5_n2_av
bsr_diff_svm_kde.5.40.5_n2_av<-BayesianSignedRank(as.matrix(diff_svm_kde.5.40.5_n2_av),-0.01,0.01)
bsr_diff_svm_kde.5.40.5_n2_av
## $winLeft
## [1] 0.1084667
##
## $winRope
## [1] 0.01506667
##
## $winRight
## [1] 0.8764667
bsr_diff_svm_kde.5.40.5_n2_av_odds.left<-bsr_diff_svm_kde.5.40.5_n2_av $winLeft/bsr_diff_svm_kde.5.40.5_n2_av $winRight
bsr_diff_svm_kde.5.40.5_n2_av_odds.left
## [1] 0.1237545
plot(rope(diff_svm_kde.5.40.5_n2_av,c(-0.01,0.01)))

diff_svm_kde.5.40.5_n3_av<-svm_dataset_av - svm_kde.5.40.5_n3_av
bsr_diff_svm_kde.5.40.5_n3_av<-BayesianSignedRank(as.matrix(diff_svm_kde.5.40.5_n3_av),-0.01,0.01)
bsr_diff_svm_kde.5.40.5_n3_av
## $winLeft
## [1] 0.06163333
##
## $winRope
## [1] 0.05376667
##
## $winRight
## [1] 0.8846
bsr_diff_svm_kde.5.40.5_n3_av_odds.left<-bsr_diff_svm_kde.5.40.5_n3_av $winLeft/bsr_diff_svm_kde.5.40.5_n3_av $winRight
bsr_diff_svm_kde.5.40.5_n3_av_odds.left
## [1] 0.06967368
plot(rope(diff_svm_kde.5.40.5_n3_av,c(-0.01,0.01)))

diff_svm_kde.5.40.5_n4_av<-svm_dataset_av - svm_kde.5.40.5_n4_av
bsr_diff_svm_kde.5.40.5_n4_av<-BayesianSignedRank(as.matrix(diff_svm_kde.5.40.5_n4_av),-0.01,0.01)
bsr_diff_svm_kde.5.40.5_n4_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.03833333
##
## $winRight
## [1] 0.9616667
bsr_diff_svm_kde.5.40.5_n4_av_odds.left<-bsr_diff_svm_kde.5.40.5_n4_av $winLeft/bsr_diff_svm_kde.5.40.5_n4_av $winRight
bsr_diff_svm_kde.5.40.5_n4_av_odds.left
## [1] 0
plot(rope(diff_svm_kde.5.40.5_n4_av,c(-0.01,0.01)))

diff_svm_kde.5.40.5_n5_av<-svm_dataset_av - svm_kde.5.40.5_n5_av
bsr_diff_svm_kde.5.40.5_n5_av<-BayesianSignedRank(as.matrix(diff_svm_kde.5.40.5_n5_av),-0.01,0.01)
bsr_diff_svm_kde.5.40.5_n5_av
## $winLeft
## [1] 0.06286667
##
## $winRope
## [1] 0.0516
##
## $winRight
## [1] 0.8855333
bsr_diff_svm_kde.5.40.5_n5_av_odds.left<-bsr_diff_svm_kde.5.40.5_n5_av $winLeft/bsr_diff_svm_kde.5.40.5_n5_av $winRight
bsr_diff_svm_kde.5.40.5_n5_av_odds.left
## [1] 0.070993
plot(rope(diff_svm_kde.5.40.5_n5_av,c(-0.01,0.01)))

######################### Neural Network
##Neural Network Results
nn1_dataset_av<-c(0.7592, 0.3779, 0.9799)
nn1_pca.5.40.5_n1_av<-c(0.2408, 0.3995, 0.8983)
nn1_pca.5.40.5_n2_av<-c(0.3060, 0.5468, 0.9523)
nn1_pca.5.40.5_n3_av<-c(0.5026, 0.3311, 0.3753)
nn1_pca.5.40.5_n4_av<-c(0.7878, 0.2647, 0.1017)
nn1_pca.5.40.5_n5_av<-c(0.7592, NA, NA)
nn1_kde.5.40.5_n1_av<-c(0.8300, 0.2059, 0.9808)
nn1_kde.5.40.5_n2_av<-c(0.7592, 0.4429, 0.9309)
nn1_kde.5.40.5_n3_av<-c(0.7893, 0.5488, 0.9570)
nn1_kde.5.40.5_n4_av<-c(0.7054, 0.3961, 0.9570)
nn1_kde.5.40.5_n5_av<-c(0.7690, 0.3801, 0.8983)
######################## ROPE PCA
diff_nn1_pca.5.40.5_n1_av<-nn1_dataset_av - nn1_pca.5.40.5_n1_av
bsr_diff_nn1_pca.5.40.5_n1_av<-BayesianSignedRank(as.matrix(diff_nn1_pca.5.40.5_n1_av),-0.01,0.01)
bsr_diff_nn1_pca.5.40.5_n1_av
## $winLeft
## [1] 0.1099667
##
## $winRope
## [1] 0.0144
##
## $winRight
## [1] 0.8756333
bsr_diff_nn1_pca.5.40.5_n1_av_odds.left<-bsr_diff_nn1_pca.5.40.5_n1_av$winLeft/bsr_diff_nn1_pca.5.40.5_n1_av $winRight
bsr_diff_nn1_pca.5.40.5_n1_av_odds.left
## [1] 0.1255853
plot(rope(diff_nn1_pca.5.40.5_n1_av,c(-0.01,0.01)))

diff_nn1_pca.5.40.5_n2_av<-nn1_dataset_av - nn1_pca.5.40.5_n2_av
bsr_diff_nn1_pca.5.40.5_n2_av<-BayesianSignedRank(as.matrix(diff_nn1_pca.5.40.5_n2_av),-0.01,0.01)
bsr_diff_nn1_pca.5.40.5_n2_av
## $winLeft
## [1] 0.26
##
## $winRope
## [1] 0.01703333
##
## $winRight
## [1] 0.7229667
bsr_diff_nn1_pca.5.40.5_n2_av_odds.left<-bsr_diff_nn1_pca.5.40.5_n2_av$winLeft/bsr_diff_nn1_pca.5.40.5_n2_av $winRight
bsr_diff_nn1_pca.5.40.5_n2_av_odds.left
## [1] 0.3596293
plot(rope(diff_nn1_pca.5.40.5_n2_av,c(-0.01,0.01)))

diff_nn1_pca.5.40.5_n3_av<-nn1_dataset_av - nn1_pca.5.40.5_n3_av
bsr_diff_nn1_pca.5.40.5_n3_av<-BayesianSignedRank(as.matrix(diff_nn1_pca.5.40.5_n3_av),-0.01,0.01)
bsr_diff_nn1_pca.5.40.5_n3_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.008833333
##
## $winRight
## [1] 0.9911667
bsr_diff_nn1_pca.5.40.5_n3_av_odds.left<-bsr_diff_nn1_pca.5.40.5_n3_av$winLeft/bsr_diff_nn1_pca.5.40.5_n3_av $winRight
bsr_diff_nn1_pca.5.40.5_n3_av_odds.left
## [1] 0
plot(rope(diff_nn1_pca.5.40.5_n3_av,c(-0.01,0.01)))

diff_nn1_pca.5.40.5_n4_av<-nn1_dataset_av - nn1_pca.5.40.5_n4_av
bsr_diff_nn1_pca.5.40.5_n4_av<-BayesianSignedRank(as.matrix(diff_nn1_pca.5.40.5_n4_av),-0.01,0.01)
bsr_diff_nn1_pca.5.40.5_n4_av
## $winLeft
## [1] 0.1070333
##
## $winRope
## [1] 0.0169
##
## $winRight
## [1] 0.8760667
bsr_diff_nn1_pca.5.40.5_n4_av_odds.left<-bsr_diff_nn1_pca.5.40.5_n4_av$winLeft/bsr_diff_nn1_pca.5.40.5_n4_av $winRight
bsr_diff_nn1_pca.5.40.5_n4_av_odds.left
## [1] 0.1221749
plot(rope(diff_nn1_pca.5.40.5_n4_av,c(-0.01,0.01)))

#diff_nn1_pca.5.40.5_n5_av<-nn1_dataset_av - nn1_pca.5.40.5_n5_av
#bsr_diff_nn1_pca.5.40.5_n5_av<-BayesianSignedRank(as.matrix(diff_nn1_pca.5.40.5_n5_av),-0.01,0.01)
#bsr_diff_nn1_pca.5.40.5_n5_av
#bsr_diff_nn1_pca.5.40.5_n5_av_odds.left<-bsr_diff_nn1_pca.5.40.5_n5_av$winLeft/bsr_diff_nn1_pca.5.40.5_n5_av $winRight
#bsr_diff_nn1_pca.5.40.5_n5_av_odds.left
#plot(rope(diff_nn1_pca.5.40.5_n5_av,c(-0.01,0.01)))
########################## ROPE KDE
diff_nn1_kde.5.40.5_n1_av<-nn1_dataset_av - nn1_kde.5.40.5_n1_av
bsr_diff_nn1_kde.5.40.5_n1_av<-BayesianSignedRank(as.matrix(diff_nn1_kde.5.40.5_n1_av),-0.01,0.01)
bsr_diff_nn1_kde.5.40.5_n1_av
## $winLeft
## [1] 0.3047
##
## $winRope
## [1] 0.1998
##
## $winRight
## [1] 0.4955
bsr_diff_nn1_kde.5.40.5_n1_av_odds.left<-bsr_diff_nn1_kde.5.40.5_n1_av$winLeft/bsr_diff_nn1_kde.5.40.5_n1_av$winRight
bsr_diff_nn1_kde.5.40.5_n1_av_odds.left
## [1] 0.6149344
plot(rope(diff_nn1_kde.5.40.5_n1_av,c(-0.01,0.01)))

diff_nn1_kde.5.40.5_n2_av<-nn1_dataset_av - nn1_kde.5.40.5_n2_av
bsr_diff_nn1_kde.5.40.5_n2_av<-BayesianSignedRank(as.matrix(diff_nn1_kde.5.40.5_n2_av),-0.01,0.01)
bsr_diff_nn1_kde.5.40.5_n2_av
## $winLeft
## [1] 0.3495333
##
## $winRope
## [1] 0.3047667
##
## $winRight
## [1] 0.3457
bsr_diff_nn1_kde.5.40.5_n2_av_odds.left<-bsr_diff_nn1_kde.5.40.5_n2_av$winLeft/bsr_diff_nn1_kde.5.40.5_n2_av$winRight
bsr_diff_nn1_kde.5.40.5_n2_av_odds.left
## [1] 1.011089
plot(rope(diff_nn1_kde.5.40.5_n2_av,c(-0.01,0.01)))

diff_nn1_kde.5.40.5_n3_av<-nn1_dataset_av - nn1_kde.5.40.5_n3_av
bsr_diff_nn1_kde.5.40.5_n3_av<-BayesianSignedRank(as.matrix(diff_nn1_kde.5.40.5_n3_av),-0.01,0.01)
bsr_diff_nn1_kde.5.40.5_n3_av
## $winLeft
## [1] 0.7973333
##
## $winRope
## [1] 0.04783333
##
## $winRight
## [1] 0.1548333
bsr_diff_nn1_kde.5.40.5_n3_av_odds.left<-bsr_diff_nn1_kde.5.40.5_n3_av$winLeft/bsr_diff_nn1_kde.5.40.5_n3_av$winRight
bsr_diff_nn1_kde.5.40.5_n3_av_odds.left
## [1] 5.149623
plot(rope(diff_nn1_kde.5.40.5_n3_av,c(-0.01,0.01)))

diff_nn1_kde.5.40.5_n4_av<-nn1_dataset_av - nn1_kde.5.40.5_n4_av
bsr_diff_nn1_kde.5.40.5_n4_av<-BayesianSignedRank(as.matrix(diff_nn1_kde.5.40.5_n4_av),-0.01,0.01)
bsr_diff_nn1_kde.5.40.5_n4_av
## $winLeft
## [1] 0.08173333
##
## $winRope
## [1] 0.1406
##
## $winRight
## [1] 0.7776667
bsr_diff_nn1_kde.5.40.5_n4_av_odds.left<-bsr_diff_nn1_kde.5.40.5_n4_av$winLeft/bsr_diff_nn1_kde.5.40.5_n4_av$winRight
bsr_diff_nn1_kde.5.40.5_n4_av_odds.left
## [1] 0.1051007
plot(rope(diff_nn1_kde.5.40.5_n4_av,c(-0.01,0.01)))

diff_nn1_kde.5.40.5_n5_av<-nn1_dataset_av - nn1_kde.5.40.5_n5_av
bsr_diff_nn1_kde.5.40.5_n5_av<-BayesianSignedRank(as.matrix(diff_nn1_kde.5.40.5_n5_av),-0.01,0.01)
bsr_diff_nn1_kde.5.40.5_n5_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.5833667
##
## $winRight
## [1] 0.4166333
bsr_diff_nn1_kde.5.40.5_n5_av_odds.left<-bsr_diff_nn1_kde.5.40.5_n5_av$winLeft/bsr_diff_nn1_kde.5.40.5_n5_av$winRight
bsr_diff_nn1_kde.5.40.5_n5_av_odds.left
## [1] 0
plot(rope(diff_nn1_kde.5.40.5_n5_av,c(-0.01,0.01)))

################################ Logistic Regression
##Logistic Rergression Results
lr_dataset_av<-c(0.8533, 0.9245, 0.9778)
lr_pca.5.40.5_n1_av<-c(0.2488, 0.5740, 0.8985)
lr_pca.5.40.5_n2_av<-c(0.3047, 0.8836, 0.9445)
lr_pca.5.40.5_n3_av<-c(0.4532, 0.6105, 0.7772)
lr_pca.5.40.5_n4_av<-c(0.7752, 0.3431, 0.1026)
lr_pca.5.40.5_n5_av<-c(0.6868, NA, NA)
lr_kde.5.40.5_n1_av<-c(0.8507, 0.8968, 0.9786)
lr_kde.5.40.5_n2_av<-c(0.8449, 0.8787, 0.9788)
lr_kde.5.40.5_n3_av<-c(0.8349, 0.7556, 0.9717)
lr_kde.5.40.5_n4_av<-c(0.8266, 0.6623, 0.9659)
lr_kde.5.40.5_n5_av<-c(0.8140, 0.5659, 0.9333)
######################## ROPE PCA
diff_lr_pca.5.40.5_n1_av<-lr_dataset_av - lr_pca.5.40.5_n1_av
bsr_diff_lr_pca.5.40.5_n1_av<-BayesianSignedRank(as.matrix(diff_lr_pca.5.40.5_n1_av),-0.01,0.01)
bsr_diff_lr_pca.5.40.5_n1_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.008933333
##
## $winRight
## [1] 0.9910667
bsr_diff_lr_pca.5.40.5_n1_av_odds.left<-bsr_diff_lr_pca.5.40.5_n1_av$winLeft/bsr_diff_lr_pca.5.40.5_n1_av $winRight
bsr_diff_lr_pca.5.40.5_n1_av_odds.left
## [1] 0
plot(rope(diff_lr_pca.5.40.5_n1_av,c(-0.01,0.01)))

diff_lr_pca.5.40.5_n2_av<-lr_dataset_av - lr_pca.5.40.5_n2_av
bsr_diff_lr_pca.5.40.5_n2_av<-BayesianSignedRank(as.matrix(diff_lr_pca.5.40.5_n2_av),-0.01,0.01)
bsr_diff_lr_pca.5.40.5_n2_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.008133333
##
## $winRight
## [1] 0.9918667
bsr_diff_lr_pca.5.40.5_n2_av_odds.left<-bsr_diff_lr_pca.5.40.5_n2_av$winLeft/bsr_diff_lr_pca.5.40.5_n2_av $winRight
bsr_diff_lr_pca.5.40.5_n2_av_odds.left
## [1] 0
plot(rope(diff_lr_pca.5.40.5_n2_av,c(-0.01,0.01)))

diff_lr_pca.5.40.5_n3_av<-lr_dataset_av - lr_pca.5.40.5_n3_av
bsr_diff_lr_pca.5.40.5_n3_av<-BayesianSignedRank(as.matrix(diff_lr_pca.5.40.5_n3_av),-0.01,0.01)
bsr_diff_lr_pca.5.40.5_n3_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.0095
##
## $winRight
## [1] 0.9905
bsr_diff_lr_pca.5.40.5_n3_av_odds.left<-bsr_diff_lr_pca.5.40.5_n3_av$winLeft/bsr_diff_lr_pca.5.40.5_n3_av $winRight
bsr_diff_lr_pca.5.40.5_n3_av_odds.left
## [1] 0
plot(rope(diff_lr_pca.5.40.5_n3_av,c(-0.01,0.01)))

diff_lr_pca.5.40.5_n4_av<-lr_dataset_av - lr_pca.5.40.5_n4_av
bsr_diff_lr_pca.5.40.5_n4_av<-BayesianSignedRank(as.matrix(diff_lr_pca.5.40.5_n4_av),-0.01,0.01)
bsr_diff_lr_pca.5.40.5_n4_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.007933333
##
## $winRight
## [1] 0.9920667
bsr_diff_lr_pca.5.40.5_n4_av_odds.left<-bsr_diff_lr_pca.5.40.5_n4_av$winLeft/bsr_diff_lr_pca.5.40.5_n4_av $winRight
bsr_diff_lr_pca.5.40.5_n4_av_odds.left
## [1] 0
plot(rope(diff_lr_pca.5.40.5_n4_av,c(-0.01,0.01)))

#diff_lr_pca.5.40.5_n5_av<-lr_dataset_av - lr_pca.5.40.5_n5_av
#bsr_diff_lr_pca.5.40.5_n5_av<-BayesianSignedRank(as.matrix(diff_lr_pca.5.40.5_n5_av),-0.01,0.01)
#bsr_diff_lr_pca.5.40.5_n5_av
#bsr_diff_lr_pca.5.40.5_n5_av_odds.left<-bsr_diff_lr_pca.5.40.5_n5_av$winLeft/bsr_diff_lr_pca.5.40.5_n5_av $winRight
#bsr_diff_lr_pca.5.40.5_n5_av_odds.left
#plot(rope(diff_lr_pca.5.40.5_n5_av,c(-0.01,0.01)))
########################## ROPE KDE
diff_lr_kde.5.40.5_n1_av<-lr_dataset_av - lr_kde.5.40.5_n1_av
bsr_diff_lr_kde.5.40.5_n1_av<-BayesianSignedRank(as.matrix(diff_lr_kde.5.40.5_n1_av),-0.01,0.01)
bsr_diff_lr_kde.5.40.5_n1_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.5800333
##
## $winRight
## [1] 0.4199667
bsr_diff_lr_kde.5.40.5_n1_av_odds.left<-bsr_diff_lr_kde.5.40.5_n1_av $winLeft/bsr_diff_lr_kde.5.40.5_n1_av $winRight
bsr_diff_lr_kde.5.40.5_n1_av_odds.left
## [1] 0
plot(rope(diff_lr_kde.5.40.5_n1_av,c(-0.01,0.01)))

diff_lr_kde.5.40.5_n2_av<-lr_dataset_av - lr_kde.5.40.5_n2_av
bsr_diff_lr_kde.5.40.5_n2_av<-BayesianSignedRank(as.matrix(diff_lr_kde.5.40.5_n2_av),-0.01,0.01)
bsr_diff_lr_kde.5.40.5_n2_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.5761333
##
## $winRight
## [1] 0.4238667
bsr_diff_lr_kde.5.40.5_n2_av_odds.left<-bsr_diff_lr_kde.5.40.5_n2_av $winLeft/bsr_diff_lr_kde.5.40.5_n2_av $winRight
bsr_diff_lr_kde.5.40.5_n2_av_odds.left
## [1] 0
plot(rope(diff_lr_kde.5.40.5_n2_av,c(-0.01,0.01)))

diff_lr_kde.5.40.5_n3_av<-lr_dataset_av - lr_kde.5.40.5_n3_av
bsr_diff_lr_kde.5.40.5_n3_av<-BayesianSignedRank(as.matrix(diff_lr_kde.5.40.5_n3_av),-0.01,0.01)
bsr_diff_lr_kde.5.40.5_n3_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.2177333
##
## $winRight
## [1] 0.7822667
bsr_diff_lr_kde.5.40.5_n3_av_odds.left<-bsr_diff_lr_kde.5.40.5_n3_av $winLeft/bsr_diff_lr_kde.5.40.5_n3_av $winRight
bsr_diff_lr_kde.5.40.5_n3_av_odds.left
## [1] 0
plot(rope(diff_lr_kde.5.40.5_n3_av,c(-0.01,0.01)))

diff_lr_kde.5.40.5_n4_av<-lr_dataset_av - lr_kde.5.40.5_n4_av
bsr_diff_lr_kde.5.40.5_n4_av<-BayesianSignedRank(as.matrix(diff_lr_kde.5.40.5_n4_av),-0.01,0.01)
bsr_diff_lr_kde.5.40.5_n4_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.03603333
##
## $winRight
## [1] 0.9639667
bsr_diff_lr_kde.5.40.5_n4_av_odds.left<-bsr_diff_lr_kde.5.40.5_n4_av $winLeft/bsr_diff_lr_kde.5.40.5_n4_av $winRight
bsr_diff_lr_kde.5.40.5_n4_av_odds.left
## [1] 0
plot(rope(diff_lr_kde.5.40.5_n4_av,c(-0.01,0.01)))

diff_lr_kde.5.40.5_n5_av<-lr_dataset_av - lr_kde.5.40.5_n5_av
bsr_diff_lr_kde.5.40.5_n5_av<-BayesianSignedRank(as.matrix(diff_lr_kde.5.40.5_n5_av),-0.01,0.01)
bsr_diff_lr_kde.5.40.5_n5_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.008333333
##
## $winRight
## [1] 0.9916667
bsr_diff_lr_kde.5.40.5_n5_av_odds.left<-bsr_diff_lr_kde.5.40.5_n5_av $winLeft/bsr_diff_lr_kde.5.40.5_n5_av $winRight
bsr_diff_lr_kde.5.40.5_n5_av_odds.left
## [1] 0
plot(rope(diff_lr_kde.5.40.5_n5_av,c(-0.01,0.01)))

#################################################### Naive Bayes
##Naive Bayes Results
nb_dataset_av<-c(0.7728, 0.8944, 0.9618)
nb_pca.5.40.5_n1_av<-c(0.2408, 0.5980, 0.8942)
nb_pca.5.40.5_n2_av<-c(0.2408, 0.8417, 0.9385)
nb_pca.5.40.5_n3_av<-c(0.7613, NA, 0.4772)
nb_pca.5.40.5_n4_av<-c(0.7592, 0.3473, 0.6929)
nb_pca.5.40.5_n5_av<-c(0.7592, NA, NA)
nb_kde.5.40.5_n1_av<-c(0.7742, 0.8485, 0.9695)
nb_kde.5.40.5_n2_av<-c(0.2408, 0.8417, 0.9385)
nb_kde.5.40.5_n3_av<-c(0.7613, NA, 0.4772)
nb_kde.5.40.5_n4_av<-c(0.7592, 0.3473, 0.6929)
nb_kde.5.40.5_n5_av<-c(0.7592, NA, NA)
######################## ROPE PCA
diff_nb_pca.5.40.5_n1_av<-nb_dataset_av - nb_pca.5.40.5_n1_av
bsr_diff_nb_pca.5.40.5_n1_av<-BayesianSignedRank(as.matrix(diff_nb_pca.5.40.5_n1_av),-0.01,0.01)
bsr_diff_nb_pca.5.40.5_n1_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.009333333
##
## $winRight
## [1] 0.9906667
bsr_diff_nb_pca.5.40.5_n1_av_odds.left<-bsr_diff_nb_pca.5.40.5_n1_av$winLeft/bsr_diff_nb_pca.5.40.5_n1_av $winRight
bsr_diff_nb_pca.5.40.5_n1_av_odds.left
## [1] 0
plot(rope(diff_nb_pca.5.40.5_n1_av,c(-0.01,0.01)))

diff_nb_pca.5.40.5_n2_av<-nb_dataset_av - nb_pca.5.40.5_n2_av
bsr_diff_nb_pca.5.40.5_n2_av<-BayesianSignedRank(as.matrix(diff_nb_pca.5.40.5_n2_av),-0.01,0.01)
bsr_diff_nb_pca.5.40.5_n2_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.0085
##
## $winRight
## [1] 0.9915
bsr_diff_nb_pca.5.40.5_n2_av_odds.left<-bsr_diff_nb_pca.5.40.5_n2_av$winLeft/bsr_diff_nb_pca.5.40.5_n2_av $winRight
bsr_diff_nb_pca.5.40.5_n2_av_odds.left
## [1] 0
plot(rope(diff_nb_pca.5.40.5_n2_av,c(-0.01,0.01)))

#diff_nb_pca.5.40.5_n3_av<-nb_dataset_av - nb_pca.5.40.5_n3_av
#bsr_diff_nb_pca.5.40.5_n3_av<-BayesianSignedRank(as.matrix(diff_nb_pca.5.40.5_n3_av),-0.01,0.01)
#bsr_diff_nb_pca.5.40.5_n3_av
#bsr_diff_nb_pca.5.40.5_n3_av_odds.left<-bsr_diff_nb_pca.5.40.5_n3_av$winLeft/bsr_diff_nb_pca.5.40.5_n3_av $winRight
#bsr_diff_nb_pca.5.40.5_n3_av_odds.left
#plot(rope(diff_nb_pca.5.40.5_n3_av,c(-0.01,0.01)))
diff_nb_pca.5.40.5_n4_av<-nb_dataset_av - nb_pca.5.40.5_n4_av
bsr_diff_nb_pca.5.40.5_n4_av<-BayesianSignedRank(as.matrix(diff_nb_pca.5.40.5_n4_av),-0.01,0.01)
bsr_diff_nb_pca.5.40.5_n4_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.03743333
##
## $winRight
## [1] 0.9625667
bsr_diff_nb_pca.5.40.5_n4_av_odds.left<-bsr_diff_nb_pca.5.40.5_n4_av$winLeft/bsr_diff_nb_pca.5.40.5_n4_av $winRight
bsr_diff_nb_pca.5.40.5_n4_av_odds.left
## [1] 0
plot(rope(diff_nb_pca.5.40.5_n4_av,c(-0.01,0.01)))

#diff_nb_pca.5.40.5_n5_av<-nb_dataset_av - nb_pca.5.40.5_n5_av
#bsr_diff_nb_pca.5.40.5_n5_av<-BayesianSignedRank(as.matrix(diff_nb_pca.5.40.5_n5_av),-0.01,0.01)
#bsr_diff_nb_pca.5.40.5_n5_av
#bsr_diff_nb_pca.5.40.5_n5_av_odds.left<-bsr_diff_nb_pca.5.40.5_n5_av$winLeft/bsr_diff_nb_pca.5.40.5_n5_av $winRight
#bsr_diff_nb_pca.5.40.5_n5_av_odds.left
#plot(rope(diff_nb_pca.5.40.5_n5_av,c(-0.01,0.01)))
########################## ROPE KDE
diff_nb_kde.5.40.5_n1_av<-nb_dataset_av - nb_kde.5.40.5_n1_av
bsr_diff_nb_kde.5.40.5_n1_av<-BayesianSignedRank(as.matrix(diff_nb_kde.5.40.5_n1_av),-0.01,0.01)
bsr_diff_nb_kde.5.40.5_n1_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.5872667
##
## $winRight
## [1] 0.4127333
bsr_diff_nb_kde.5.40.5_n1_av_odds.left<-bsr_diff_nb_kde.5.40.5_n1_av $winLeft/bsr_diff_nb_kde.5.40.5_n1_av $winRight
bsr_diff_nb_kde.5.40.5_n1_av_odds.left
## [1] 0
plot(rope(diff_nb_kde.5.40.5_n1_av,c(-0.01,0.01)))

diff_nb_kde.5.40.5_n2_av<-nb_dataset_av - nb_kde.5.40.5_n2_av
bsr_diff_nb_kde.5.40.5_n2_av<-BayesianSignedRank(as.matrix(diff_nb_kde.5.40.5_n2_av),-0.01,0.01)
bsr_diff_nb_kde.5.40.5_n2_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.0093
##
## $winRight
## [1] 0.9907
bsr_diff_nb_kde.5.40.5_n2_av_odds.left<-bsr_diff_nb_kde.5.40.5_n2_av $winLeft/bsr_diff_nb_kde.5.40.5_n2_av $winRight
bsr_diff_nb_kde.5.40.5_n2_av_odds.left
## [1] 0
plot(rope(diff_nb_kde.5.40.5_n2_av,c(-0.01,0.01)))

#diff_nb_kde.5.40.5_n3_av<-nb_dataset_av - nb_kde.5.40.5_n3_av
#bsr_diff_nb_kde.5.40.5_n3_av<-BayesianSignedRank(as.matrix(diff_nb_kde.5.40.5_n3_av),-0.01,0.01)
#bsr_diff_nb_kde.5.40.5_n3_av
#bsr_diff_nb_kde.5.40.5_n3_av_odds.left<-bsr_diff_nb_kde.5.40.5_n3_av $winLeft/bsr_diff_nb_kde.5.40.5_n3_av $winRight
#bsr_diff_nb_kde.5.40.5_n3_av_odds.left
#plot(rope(diff_nb_kde.5.40.5_n3_av,c(-0.01,0.01)))
diff_nb_kde.5.40.5_n4_av<-nb_dataset_av - nb_kde.5.40.5_n4_av
bsr_diff_nb_kde.5.40.5_n4_av<-BayesianSignedRank(as.matrix(diff_nb_kde.5.40.5_n4_av),-0.01,0.01)
bsr_diff_nb_kde.5.40.5_n4_av
## $winLeft
## [1] 0
##
## $winRope
## [1] 0.03606667
##
## $winRight
## [1] 0.9639333
bsr_diff_nb_kde.5.40.5_n4_av_odds.left<-bsr_diff_nb_kde.5.40.5_n4_av $winLeft/bsr_diff_nb_kde.5.40.5_n4_av $winRight
bsr_diff_nb_kde.5.40.5_n4_av_odds.left
## [1] 0
plot(rope(diff_nb_kde.5.40.5_n4_av,c(-0.01,0.01)))

#diff_nb_kde.5.40.5_n5_av<-nb_dataset_av - nb_kde.5.40.5_n5_av
#bsr_diff_nb_kde.5.40.5_n5_av<-BayesianSignedRank(as.matrix(diff_nb_kde.5.40.5_n5_av),-0.01,0.01)
#bsr_diff_nb_kde.5.40.5_n5_av
#bsr_diff_nb_kde.5.40.5_n5_av_odds.left<-bsr_diff_nb_kde.5.40.5_n5_av $winLeft/bsr_diff_nb_kde.5.40.5_n5_av $winRight
#bsr_diff_nb_kde.5.40.5_n5_av_odds.left
#plot(rope(diff_nb_kde.5.40.5_n5_av,c(-0.01,0.01)))